Microscopic origin of the adiabatic change of magnetization in molecular magnets 
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A microscopic model of the molecular magnet V15 is used to study mechanisms for the adiabatic 
change of the magnetization in time-dependent magnetic fields. Effects of the Dzyaloshinskii-Moriya 
interaction, the most plausible source for the energy-level repulsions that lead to adiabatic changes 
of the magnetization, are studied in detail. We find that the energy-level repulsions that result from 
this interaction exhibit a strong dependence on the direction of the applied field. We also discuss 
the role of magnetic anisotropy in the molecule Mni2-acetate. 

PACS numbers: 75.10Jm, 75.50.Xx; 75.45.+J; 75.50.Ee 



I. INTRODUCTION 



Recently, magnetic molecules such as Mni2 or V15 have attracted a lot of interest. These nanomagnets are often used 
to study explicit real-time quantum dynamics, e.g . tunne ling of the magnetization and quantum (de)coherence P-IUS 
IS, IM IH (H 1^ E3, US, EM EH E3, EM Asa result of the very weak intermolecular interactions, 

experiments can directly probe the magnetization dynamics of the individual molecules. In particular, the adiabatic 
change of the magnetization at low-temperature is governed by the discrete energy-level structure [23|, |2J, |25|, |2|| |23 • 

The adiabatic change of the magnetization requires some interactions that yield energy level repulsions, i.e. inter- 
actions that do not commute with the magnetization. The Dzyaloshinskii-Moriya (DM) interaction is the most likely 
candidate for such an interaction [Til l28l |2{^| . In the case of anisotropic high-spin molecules such as Mni2 and Fes, 
simplified anisotropic single-spin models for a specific spin multiplet can approximately reproduce the gaps of the level 
repulsions. However, the case of V15 is more complicated because V15 has half-odd-integer spin and the time-reversal 
symmetry enforces at least a two-fold degeneracy of the energy levels at zero field. 

It has been pointed out that the DM interaction is accompanied by a higher-order correction term that restores the 
SU(2) symmetry [30l l3ll I32I l33| ■ In this paper, we focus on the effects of the DM interaction and leave the inclusion 
of the higher-order term for future study. As another source of level repulsion, we might consider the hy per fi ne 
interaction with the nuclear spin. The effects of the hyperfme interaction have already been discussed in Ref. 

As the DM interaction has a vector character and is anisotropic, the dynamics of the magnetization is expected to 
depend on the direction of the magnetic field. First we study the characteristic properties of the DM interaction for 
a simplified model of V15, namely three spins on a triangle. Then we confirm the properties found in the three-spin 
model by full diagonalization of the 15-spin model of V15. In order to bridge the energy scales involved (e.g. from 
800K, a typical energy scale for the interaction between individual magnetic ions, to about 10 _2 K, a typical energy 
scale for energy- level splittings), the calculation of the energy levels of the many-spin Hamiltonian has to be very 
accurate. We have tested various standard algorithms to compute the low-lying states. For systems that are too 
large to be solved by full exact diagonalization (such as the 15-spin V15 model), we use the Lanczos method with 
full orthogonalization (LFO), a Chebyshev polynomial projector (CPP) method, and a power method with additional 
subspace diagonalization. These algorithms can solve the rather large eigenvalue problems with sufficient accuracy. 
The consistency of the data obtained by different methods gives extra confidence in the numerical results. 

The magnetic properties of molecules such as Mni2 are often studied by considering a simplified model for the 
magnetic energy levels for a specific spin multiplet, e.g. S = 10. However, for these and other, similar, magnetic 
molecules that consist of several magnetic moments (in the case of Mni2i eight Mn 3+ (S f = 2) and four Mn 4+ (S f = 
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3/2)), the reduction of the many-body Hamiltonian to an effective Hamiltonian for a specific spin multiplet is non- 
trivial. Magnetic anisotropy, a result of the geometrical arrangement of the magnetic ions within a molecule of 
low symmetry, mixes states of different total spin and enforces a treatment of the full Hilbert space of the system. 
For Mni2, the dominant contribution to the magnetic mixing due to spin-orbit interactions is also given by the 
DM interaction |35l |36(. In principle, this type of interaction can change energy-level crossings into energy-level 
repulsions. The presence of the latter is essential to explain the adiabatic changes of the magnetization at the 
resonant fields |23l l24l l25l l26t l27j . Thus, a minimal magnetic model Hamiltonian should contain (str ong ) Heisenberg 
interactions, DM interactions and a coupling to the applied magnetic field |llj, |28|, l29l 1371 l38l I39L |4C1 l4ll |42j, |43| . 
Experiments on Mni2 suggest that the energy gaps related to the transition from a state with magnetization M w — 10 
to a state with magnetization M < 4 are of the order of 10 _9 K |44|. Such gaps are too small to detect with standard 
precision (13-14 digits) calculations and therefore, in this paper we only present the global energy level diagram 
obtained from microscopic model calculations. 

The paper is organized as follows. In Sec. ^ we analyze a reduced three-spin model with C3 symmetry for the Vj.5 
molecule. Results for the energy level schemes for a 15-spin model of the V15 are presented in Sec. lIIII In Sec. lIVI we 
discuss the effects of anisotropic terms and the cases with less symmetry. In Sec. El we report results for Mni2. In 
Sec. I VII we give our conclusions. In Appendix A and B, we briefly discuss some analytical solutions and the numerical 
algorithms that we use to compute the energy levels, respectively. 



II. TRIANGLE-MODEL WITH C 3 SYMMETRY 



As a simplified model for the V15 molecule we consider a system of three spins on a triangle with C3 symmetry 
We choose the z-axis to lie along the axis of C3 symmetry. The Hamiltonian is given by 



J2 JiA+iSi ■ Si+i + J2 D • t S * x S j] - h ' fe S * J > 

j=l (id) V i J 



n = - > , j m+ iS< ■ s i+1 + ^ d w • [Si x s,] - h • ( > j Si ) , (1) 

(id) 

where = J2.3 = Jz,i = J, denotes the exchange interaction and h represents the applied magnetic field. In general 
we can choose any direction of the DM vector Djj unless there is some additional symmetry. In the present case, 
because of the C3 symmetry, the ^-component of the DM vectors must all be equal, i.e. 

£1,2= DZ i3 =I% tl =D t , (2) 

and the x and y-components of the DM vector, Dfj, D?j, must obey the relation 



u i,z \ _ / — 5 ~~T I / x \ \ 3,1 





(3) 



where D x = Df 2 and D y = D\ 2 . 

In Appendix A we give the analytic expressions for the eigenvalues and eigenvectors of model JIJ in the case that the 
applied field is parallel to the z-axis and the conditions (0) and © hold. Then, the Hamiltonian 0} is block diagonal, 
the matrix containing 4 blocks of 2 x 2 matrices. The (8-dimensional) Hilbert space separates into 4 two-dimensional 
spaces 



{|3/2,3/2>»} , {|3/2,l/2),|a>} , {|3/2, -1/2), |6)} , {|3/2, -3/2), |&», (4) 

where the state \S, M) denotes an eigenstate of the Heisenberg model with total spin S and magnetization M. The 
expressions for the orthonormal states \a) and \b) are given in appendix A. The state \a) (\b)) denotes the state \a) 
(| 6)) with all spins reversed. From Eq. J3J it follows that there is no mixing among the four levels with S = ±1/2, 
and therefore they simply cross each other. From the analytical solution it is also easy to see that there are no energy 
level repulsions at h = 0. Furthermore one can show analytically that it is impossible to change the magnetization 
from —3/2 to +3/2 by (adiabatically) reversing the external field along the z-axis. 

In the left panel of Fig. [l]we show the energy levels as a function of the strength h of the applied magnetic field 
for the case D x = D y = D z . The field is aligned along the z-axis. It is evident that there is no level mixing at ft- = 0. 
It is important to note that the small energy difference between the first excited state and the degenerate ground 
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FIG. 1: Energy levels of Hamiltonian Q for J = -2.5K, D x = D y = D z = 0.25K. Left: Applied magnetic field h parallel to 
the z-axis. Right: Applied magnetic field h along the :r-axis. 
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FIG. 2: Left: Energy levels of Hamiltonian Jy for J = — 2.5K, D x = D y = D z = 0.25K and the applied magnetic field 
h = h(l, 0, l)/\/2 tilted by 45° with respect to the z-axis. Right: Detailed view of the /i-dependence of the four lowest energy 
levels. 



states at h = is not a tunneling gap. Therefore the magnetization does not change if the field h changes sign. On 
the other hand, if we apply the magnetization along the x-direction (i.e. perpendicular to the symmetry axis of the 
triangle), then two of the states become degenerate. We find two nearly degenerate avoided level crossings, as shown 
in the right panel of Fig.^ As an intermediate case, in Fig. [3 we show the energy level diagram for the case that the 
magnetic field is tilted by 45° with respect to the z-axis (h = ft.(l,0, l)/\/2)- Then, there is a simple crossing at zero 
field and an avoided level crossing between the levels of M — ±1/2 at a finite magnetic field, fndeed, a closer look 
at the level diagram (see right panel in Fig. [2Jl reveals that the minimum energy difference between the two pairs of 
levels does not occur at zero field but at h ~ 0.35T. This implies that the Landau-Zener-Stiickelberg transition from 
the 1 1/2 , —1/2) to the 1 1/2, 1/2) level does not take place at h = but at h « 0.35T. In conclusion, the position and 
the energy splitting of the avoided level crossing that is responsible for the adiabatic change depend on the direction 
of the field. 

The numerical results discussed above have been obtained for D x = D y = D z . In Ref. (2^ the DM vector is taken 
parallel to the ?/-axis at all the bonds and the field is applied along the z-axis. For the present model, this case 
corresponds to the case with only D z (D x = D y = 0) and a field applied in the x-direction. In this case, the gap 
opens symmetrically with respect to the applied field |22|. However, Ref. [23] did not address the dependence on the 
direction of the magnetic field. In the next section we repeat the analysis of this dependence for a 15-spin model of 
the Vis molecule. 
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FIG. 3: Schematic diagram of the magnetic interactions in model £Q| of the V15 molecule. 



III. 15-SPIN MODEL FOR THE VANADIUM COMPLEX V i5 

A. Spin interactions in V15 

In Fig. |3 we show the schematic diagram of the dominant magnetic (Heisenberg) interactions in the 15-spin model 
of the V15 molecule (K 6 [V(5'As6042(H20)]-8ri20). The magnetic structure consists of two hexagons with six S = 1/2 
spins each, enclosing a triangle with three S — 1/2 spins. All dominant Heisenberg interactions are antiferromagnetic. 
The dimension of the Hilbert space of this model is 2 15 = 32768. The minimal Hamiltonian is given by expression 
Q with 15 instead of 3 spins |22l l28l 12^. |40| . The Heisenberg interactions Ji.j in Eq. Q between the vanadium 
atoms are defined according to Fig. [3J For simplicity, we assume that D^.j = except for bonds for which the 
Heisenberg exchange constant is the largest (i.e. equal to J ) Hjl^cj. Rotations about 27r/3 and 47r/3 around the 
axis perpendicular to and passing through the center of the hexagons leave the V15 complex invariant. This enforces 
the constraints @ and 10) on the values of Djj [2!J, |43j ■ 

B. Energy level diagrams 



First we discuss our results for the lowest energy levels for various sets of estimates for the model parameters of the 
15-spin model of V15 HEIUES- For thc model parameters given in Ref. g3> J = -800K, J X =J' = -54.4K, and 
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FIG. 4: Left: The eight lowest energy levels of V15 model (0 with model parameters taken from Ref. ^3 as a function of the 
applied magnetic field h parallel to the z-axis. The values of the DM vector are given in the text. Right: Detailed view of the 
four lowest energy levels at ft w 0. 



TABLE I: The eight lowest eigenvalues Ei and total spin Si of the corresponding eigenstates of the V15 model Q with model 
parameters taken from Ref. |40|| for two values of the external applied field h parallel to the z-axis. The values of the DM vectors 
are given in the text. The distance between E t and the exact eigenvalue closest to Ei is Ai = {tpi\(H — i?j) 2 |( / 3 i ) 1/ ' 2 < xlCP 9 
for i = 1, . . . , 7. 
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J2 = J" = — 160K, J3 = J4 = J5 = Je = and in the absence of DM interactions, we find for the energy gap between 
the ground state and the first excited state at h — a value of 4.12478K. This reproduces the result of Ref. 01 an d 
is also in reasonable agreement with the experimental value of approximately 3.7K |44j . 

Following Ref. [43j we now take for the DM interaction parameters D'f 2 — D\ 2 — Df 2 — 40K, which is ap- 
proximately 5% of the largest Heisenberg coupling. Using the rotational symmetry of the hexagon it follows from 
Eq.@ that Df j4 = 14.641K, D| 4 = -54.641K, D| 4 = 40K and Df 6 = -54.641K, D y 5 6 = 14.641K, Df 6 = 40K. 
If the two hexagons are not equivalent we cannot use the symmetry to reduce the number of free parameters (see 
below). However, for simplicity, we may assume that the (2, y)-positions of the spins on the lower hexagons dif- 
fer from those on the upper hexagon by a rotation about 7r/3. This yields for the remaining model parameters 
D x W m = -14.641K, Z^ 0jll = 54.641K, Df 0jll = 40K, Df 2A3 = -40K, D\ 2 ^ = -40K, Df 2>13 = 40K, and 
Df 4 15 = 54.641K, D\ A 15 = — 14.641K, Df 4 15 = 40K. For this choice of model parameters, the eight lowest en- 
ergies of V15 model Q for two values of the applied magnetic field (h = and h = AT) along the z-axis can be 
found in Table D From Table [Q we see that for zero field, the DM interaction splits the doubly-degenerate doublets 
of S = 1/2 states into two doublets of S = 1/2 states. The difference in energy between the doubly-degenerate, first 
excited states and the two-fold degenerate ground states has a value of 0.0085K. This value is much smaller than the 
experimental estimate of 0.05K [22j, but of the same order of magnitude as the value cited in Ref. 29]. The next four 
higher levels are 5 = 3/2 states. The energy-level splitting between the 5 = 3/2 and 5=1/2 states is 4. IK. The 
value of this splitting weakly depends on the value of the DM interaction (results not shown). 

Following Ref. [if, we take J = -800K, J x = J' = -225K, J 2 = J" = -350K, and J 3 = J 4 = J s = J 6 = 0. In the 
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absence of DM interactions, we find that the energy gap between the four-fold degenerate ground state and the first 
excited state is 3.61K, in full agreement with the result of Ref. [2^. Note that this value of the gap is fairly close to 
the experimental value of 3.7K 44]. Taking for the DM interactions Df 2 = Df 415 = 25K, D| 4 = Df 6 = Df n = 
Df 213 = — 12.5K, Z?f 4 = — D\q — —D\ Q11 — Df 213 = — 21.5K, our calculation for the splitting between the two 
doubly-degenerate S — 1/2 levels yields 0.0037K, about a factor of two larger than the value cited in Ref. [2^. For 
the energy splitting between the S = 1/2 and S — 3/2 levels we obtain 3.616K instead of the value 3.618K given in 
Ref. [2j| . These differences seem to suggest that a perturbation approach [2£j for the DM interaction may not be 
sufficiently accurate for quantitative purposes 45] . 

The most advanced estimation of the model parameters is given in Ref. [l8J. Taking J = — 809K, J' = — 120K, 
J" = 120K, Ji = -30K, J 2 = -122K, J 3 = -3K, J 4 = —UK, J 5 = -3K, J 6 = -2K (see Table I in Ref. 0) yields an 
energy gap of 4.915K, in agreement with Ref. 18]. Adding the DM interaction with parameters Df 2 = Df 4 15 = 25K, 
D\ 4 = Df 6 = Df Q11 — Df 213 — — 12.5K, D 34 = — _D| 6 = — D\ n = D\ 2 13 = — 21.5K, a level repulsion appears 
at h w 3.6T where the S = 1/2 and S — 3/2 states mix (results not shown). Although the qualitative features of 
the energy-level diagram for this set of model parameters also agree with what one would expect on the basis of 
experiments, the field at which the states 1 1/2, 1/ 2} and |3/2, 3/2) cross, h 3.6T, does not compare well to the 
experimental estimate h 2.8T. 

In Fig. H we show the results for J = -800K, J x = J 1 = -225K, and J 2 = J" = -350K 29] and the same DM 
interaction parameters as those used to obtain the results of Table [I] For the energy gap at zero field, we find 
3.7K (instead of 4. IK for the J's of Ref. H3]), in good agreement with the experimental estimate of 3.7K [44|. The 
transition between the states |l/2,l/2) and |3/2, 3/2) takes place at h w 2.8T (instead of h « 3.0T for the J's of 
Ref. ^3)' a l so m good agreement with the experimental value 2.8T. 

The results for the zero-field energy gap suggest that there are many different sets of model parameters that 
approximately reproduce the experimental gap between the singlet and triplet states. However, as we show below, the 
energy gap at zero field not necessarily corresponds to the gap of a level repulsion that is required for the adiabatic 
change of the magnetization. 

We now focus on the relation between the direction of the magnetic field and the structure of the energy levels. 
As shown in the right panel of Fig.0] when the field is parallel to the symmetry axis the energy levels of M = ±1/2 
simply cross, just as in the case of the left panel of Fig.^ In Fig. we present results for the cases where the angle 
between the applied field and the z-axis is 45° and 90° , respectively. Clearly, we find the same type of dependence of 
the energy levels on the angle of the field as in the case of the three-spin model (see right panel of Fig. [21 and right 
panel of Fig. ^ respectively). Exactly the same qualitative features are obtained for the other sets of parameters 
discussed above (results not shown). 

Up to now, we used DM vectors that satisfy the rotational symmetry of a hexagon, and we took the same DM 
vectors for the other hexagon for simplicity |29j| . However, if there is some symmetry which connects the upper and 
lower hexagon, we have a relation between the DM vectors on both hexagons. In concert with the relations between 
the exchange couplings, let us assume that the upper and lower hexagon are related to each other by a 180° rotation 
around a vector that passes through V-atom number 7 (see Fig. |3J) and the middle of the line connecting the two 
other V-atoms of the triangle. This symmetry operation X transforms the sites {(1,2), (3,4), (5,6)} into {(14,15), 
(10,11), (12,13)}. If we place the y-axis along the line through V-atom number 7 and through the middle of the line 
connecting V-atoms 8 and 9, and if we take Di 2 = (D x , D y , D z ) as the reference DM vector, the other DM vectors 
are given by D 3i4 = R 2 T) ly2 , D 5i6 = i?Di i2 , D 14;ls = XD li2) D1041 = RXT) ly2 , and D 12 ,i3 = i? 2 lD lj2 . Here R 
denotes a rotation of the hexagons around 27r/3 in the plane of the hexagons. In Figs. El and [7| we show the energy 
level diagrams for the set of parameters taken from Ref. [401 ] , with the additional constraint imposed by the symmetry. 
Qualitatively the results are similar to those obtained in the absence of a symmetry that relates the upper and lower 
hexagons (see Figs. 0] and 0. 

Summarizing: As in the case of the three-spin model, all our results for the 15-spin V15 model clearly demonstrate 
that the mixing of levels strongly depends on the direction of the magnetic field. It seems therefore that this dependence 
is a generic feature of the DM interaction. 

IV. EFFECTS OF LOWER SYMMETRY 

We now discuss the effects of distortion of the triangle and anisotropic exchange interactions in the triangle-model 
of V15 on the energy level diagram for /i«0. When the triangle is distorted ( / ^2,3 / J3.1 / Jx,i) the degeneracy 
of the two doublets at h = is lifted, even if D x = D y = D z = 0. In Figs. 00 and I^Jwe show data for J\, 2 — — 2.5K, 
J 2 , 3 = -2.0K and J 3A = -3.0K and D x = D y = D z = 0.25K for the applied magnetic field h = h(l, 0, 1)/V2 tilted 
by 45° with respect to the z-axis, and for field directions parallel to the z-axis and along the £-axis. Unless the field 
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FIG. 5: Same as the right panel in Fig.[I]except that the applied magnetic field h = h(l, 0, l)/v2 is tilted by 45° with respect 
to the z-axis (left) and h is along the z-axis (right). 
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FIG. 6: Left: The eight lowest energy levels of V15 model Q with model parameters taken from Ref. |4oT | assuming that the 
DM vectors of the upper and lower hexagon are related to each other by a 180° rotation around a vector that passes through 
V-atom number 7 (see Fig. |3J of the triangle and the middle of the line connecting the two other V-atoms of the triangle. The 
applied magnetic field h is parallel to the z-axis. Right: Detailed view of the /i-dependence of the four lowest energy levels. 



is parallel to the z-axis, the energy level diagrams are qualitatively similar to the ones of the undistorted triangle (see 
Figs- Gland EJ), that is, crossings at h = and level repulsions close to h = but not at ft, = 0. However, in the case 
of an undistorted triangle with anisotropic, antiferromagnetic exchange interactions (i.e. different for the x, y, and 
z-components of the spins) and D x = D y = D z = 0.25K, the energy level diagrams are qualitatively similar to the 
ones of the undistorted triangle (results not shown) . 



V. MANGANESE COMPLEX: Mn 12 

The four inner Mn+ 4 ions in the Mni 2 molecule (Mni 2 (CH 3 COO)i 6 (H 2 0)40i2-2CH3COOH-4H 2 0), have spin 
S = 3/2. The other eight Mn +3 ions have spin S — 2. The dimension of the Hilbert space of this system is 
4 4 x 5 8 = 10 s . If the total magnetization is a conserved quantity, it can be used to block-diagonalize the Hamiltonian, 
allowing a numerical study of models of this size [4l],|4(|. However, to study the adiabatic change of magnetization, 
we have to take into account all the states, and the dimension of the matrix becomes prohibitively large. Thus we 
need to simplify the model in order to reduce the dimension. A drastic reduction of the number of spin states can 
be achieved by assuming that the strong antiferromagnetic Hcisenberg interaction J' between an inner ion and its 
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FIG. 7: Same as the right panel in Fig. HJ except that the applied magnetic field h = h(l, 0, l)/v2 is tilted by 45° with respect 
to the z-axis (left) and h is along the x-axis (right). 
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FIG. 8: Energy level diagram of the Hamiltonian Q of a distorted triangle with model parameters Ji,2 = — 2.5K, J2,3 = —2. OK, 
J3,i = —3. OK, and D x — D y — D z = 25K. The applied magnetic field h = h(l, 0, l)/v2 is tilted by 45° with respect to the 
z-axis. 



outer neighbor allows the replacement of the magnetic moment of an inner ion by an effective S — 1/2 moment. The 
schematic diagram of this simplified (but still complicated) model is shown in the left panel of Fig. ^| The 
dimension of the Hilbert space of this model is 2 4 x 5 4 = 10 4 . In the following we study this simplified model. 
The Hamiltonian for the magnetic interactions of the simplified Mni2 model can be written as [3]} 



4 2 4 8 

n = -J(E S ^) - J ' E s. i .s,-K z ^(^ +4 ) 2 + Y, r),, s, s; X> s <- ( 5 ) 

»=1 B i=l (i,j)e B i=l 

where the index l<i<4(5<i<8) refers to S = 1/2 (S = 2) spins and B denotes the set of pairs B = 
{(1, 5), (1, 8), (2, 5), (2, 6), (3, 6), (3, 7), (4, 6), (4, 8)}. The first two terms describe the isotropic Heisenberg exchange 
between the spins. The third term (K z ) describes the single-ion easy-axis anisotropy of S = 2 spins. The fourth term 
represents the anti-symmetric DM interaction in Mn.12. The vector Djj determines the DM interaction between the 
i-th S — 1/2 spin and the j-th. S — 2 spin. The last term describes the interaction of the spins with the external 
magnetic field h. Note that the factor gfis is absorbed in our definition of h. Model (J5J reproduces experimental data, 
such as the splitting of the neutron scattering peaks, results of EPR measurements and the temperature dependence 
of the magnetic susceptibility [37j ■ 
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FIG. 9: Same as Fig. [SJexcept that the applied magnetic field h is parallel to the z-axis (left) and along the :r-axis (right). 
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FIG. 10: Left: Schematic diagram of the magnetic interactions of the simplified model © of the Mni2 molecule Wfy . Black 



circles: S = 1/2; open circles: 5 = 2. Also shown are the DM vectors (for i < j and Di.j 



D 7 



Right: The 21 lowest 



energy levels of the Mni2 model as a function of the applied magnetic field h. Solid lines: Eigenstates with S « 10; dashed 
lines: Eigenstates with S ~ 9. The applied magnetic field h is parallel to the z-axis. 



The first three terms in Hamiltonian J^J conserve the z-component of the total spin M z = J2i=i The DM inter- 
action, on the other hand, mixes states with different total spin. Hence, the DM interaction can change level crossings 
into level repulsions. Therefore, the presence of the DM interaction may be sufficient to explain the experimentally 
observed adiabatic changes of the magnetization. 

The four- fold rotational-reflection symmetry (S4) of the Mni2 molecule imposes some relations between the DM- 
vectors. There are only three independent DM-parameters 37]: D x = Df s , D y = Df g , and D z = D\ 8 , as indicated 
in Fig. The parameters of model (J5J have been estimated by comparing experimental and theoretical data, fn 
this paper we will use the parameter set B from Ref. j33, 0: J = 23.8K, J' = 79. 2K, K z = 5.72K, D x = 22K, 
D y = 0, and D z — 10K. Although the amount of available data is not sufficient to fix all these parameters accurately, 
we expect that the general trends in the energy-level diagram will not change drastically if these parameters change 
relatively little. 

In Table ITTl we present results for the energy and total spin of the 21 lowest states for h = and h = 5T (in these 
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calculations h is parallel to the z-axis). The numerical results obtained by full exact diagonalization (LAPACK), the 
Lanczos method with full orthogonalization (LFO) (see appendix B) and the Chebyshev-polynomial-based projector 
method (CPP) (see appendix B) are the same to working precision (about 13 digits). Clearly there are states with 
total spin 8, 9 and 10 within these 21 lowest eigenstates. Although the total magnetization is not a good quantum 
number, we can label the various eigenstates by their (calculated) magnetization. 
The S — 10 single-spin model for Mni2 



H = -D{S z f-hS z , (6) 

where D denotes the uniaxial anisotropy, is often used as a starting point to interpret experimental results [|| 0, 0, 
IT3.llHl4fl| . The energy levels of this model exhibit crossings at the resonant fields h = ±Dn for n = —10, . . . , 10, in 
qualitative agreement with our numerical results (shown in the right panel of Fie. 110(1 for the microscopic model 0. 
A fit of the first eight level crossings of model © to the data of Fig. 1101 yields D ss 0.74K, in good agreement with 
experiments 0, 0] • The Hamiltonian of the single-spin model © commutes with the magnetization S z and therefore 
its energy diagram only displays level crossings, no level repulsions. Adding anisotropy of the form C(S+ + St) 
to mod el (B p) changes the estimated value of D and leads to level repulsions when the magnetization changes by 4 

It is also of interest to compare the level splitting at h = obtained b y lo west-order degenerate perturbation theory 
of model © with fourth-order anisotropy of the form St + St |48l l5ll l52l ] 



AF -<>2Df C Y l/2 1 (5 + m)! (7) 

AE l+1J -S2D^—j ((m/2 _ 1)!)2 (5 _ m)! - ( ? ) 

for m even (AEj+i^ = for m odd) with the result of the microscopic model calculation based on model J5J. In 
Eq. (J7J I denotes the perturbed eigenstates in increasing order of energy and m is the magnetic quantum number of 
the unperturbed states |4Sl]. Using the values D « 0.69K, C/D 5.7 x 10~ 5 obtained by fitting the single-spin model 
to experimental data |48j, the energy gap for m = 6 is given by 



AE 13 . 12 w 0.00022X. (8) 

Taking into account that (because of the presence of S w 9 states) the corresponding S ~ 10 levels of model (JjjJ are the 
14-th and 15-th lowest energy level, TableHIlshows that for h = 0, AE 13il2 = 0.00042K. In view of the uncertainties 
on the estimates of the various model parameters, the difference of only a factor of two is remarkably small. From this 
comparison, we may conclude that the DM interaction leads to energy gaps that are of the same order of magnitude 
as the gaps due to the fourth-order terms + St in the single-spin model. 

In the right panel of Fig. ^| we show the results for the lowest 21 energy levels of the Mni2 model as a function of 
the applied magnetic field as obtained by LFO. The applied magnetic field is parallel to the z-axis. In Fig. 1101 solid 
(dashed) lines represent eigenstates with S « 10 (9) (within an error of about 10%). Also eigenstates with S « 8 
appear for h > 4 but these are not shown for clarity. For the Mni2 model, the DM induced energy splittings between 
the S s« 10, M w — 10 state and other states are less than 10 _6 K. Adding an extra transverse field by tilting the 
h-field by 5° with respect to the z-axis does not change this conclusion (results not shown). 

Finally, we added to model © the next-to-lowest order relativistic correction to the local anisotropy that is com- 
patible with the symmetry of the square [37j 



Hi = K x [{Sff + {Slf + (Sf) 2 + (SI) 2 ]. (9) 

Although we took a perhaps unrealistically large value of K\ [K\ = K z /2), we were unable to detect energy level 
repulsions up to the M w —10, M w 3 transition (results not shown). On the other hand in experiments [aTl?! Il2|. 
adiabatic changes of the magnetization have been observed at h w 3.4T (M z sa — 10 — > M, « 4) and h « 3.9T 
(M z w — 10 — * M z 3) and the magnitude of the energy splittings is of the order of 10 nK |44|. The precision of the 
present calculations is about 10 _6 K. Thus, it is consistent that within the (very high) resolution in the h-ficld and 
the 13-digit precision of the calculation, no information about the gap could be extracted. The algorithms developed 
for the work presented in this paper can be used for 33-digit calculations without modification and we leave the 
calculation of the splittings for future work. 
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TABLE II: The 21 lowest eigenvalues Ei and total spin Si of the corresponding eigenstates of the Mni2 model for two 
values of the external applied field h along the z-axis. The distance between Ei and the exact eigenvalue closest to Ei is 
At = {<pi\(H - Ei) 2 \(fi) 1/2 < 1(T 10 for i = 1, ... ,7. Note that for h = the levels 12,13,18,19, and 20 belong to the S « 9 
subspace and not to the 5* ~ 10 subspace. 



i 


Ei(h = 0) 


Si(h = 0) 


Ei(h = 5T) 


5, (h = ! 





-812.771882673675 


9.72 


-878.203468556749 


9.77 


1 


-812.771882673460 


9.72 


-857.042137859145 


9.78 


2 


-798.326618260922 


9.72 


-837.727273846391 


9.76 


3 


-798.326618261122 


9.72 


-820.218084451590 


9.73 


4 


-785.677644659194 


9.70 


-816.530388063144 


8.82 


5 


-785.677644658983 


9.70 


-804.449056631056 


9.69 


6 


-774.774953284432 


9.68 


-804.242954124067 


8.77 


7 


-774.774953284294 


9.68 


-798.519706376890 


8.82 


8 


-765.549187817101 


9.65 


-790.336285910398 


9.65 


9 


-765.549187333902 


9.65 


-788.717912654407 


8.78 


10 


-757.919915510036 


9.61 


-782.037906282298 


8.82 


11 


-757.919915509970 


9.61 


-777.789121186874 


9.60 


12 


-757.673722613912 


8.77 


-774.222534072884 


8.79 


13 


-757.673722613981 


8.77 


-773.614730624955 


8.80 


14 


-751.806498916496 


9.57 


-766.993852023410 


8.80 


15 


-751.806072514140 


9.57 


-766.720405060281 


9.55 


16 


-747.135398548595 


9.54 


-760.807893152451 


8.79 


17 


-747.135398548602 


9.54 


-760.482768227423 


8.12 


18 


-746.357522623039 


8.77 


-757.060762637193 


9.51 


19 


-746.357522623082 


8.77 


-754.700878489864 


8.59 


20 


-745.778951523327 


8.70 


-753.310517350023 


8.78 



VI. DISCUSSION 



We have studied the dependence of the energy level diagrams, with level repulsions due to the Dzyaloshinskii- 
Moriya interaction, on the direction of the applied magnetic field. We found that the dependence on the direction 
of the magnetic field seems to be generic, at least if the system has C3 symmetry. Our numerical data suggest that 
the three-spin model reproduces the main features of the 15-spin model of V15. The presence of the Dzyaloshinskii- 
Moriya interaction allows for adiabatic changes of the magnetization but, according to our calculations, the value 
of the resonant field for the 1 1/2, — 1/2) to 1 1/2, 1/2) transition changes with the direction of the magnetic field. 
The Dzyaloshinskii-Moriya interaction not only lifts the degeneracy but, depending on the direction of the field with 
respect to the symmetry axis, also shifts the resonant point away from h = 0. 

The butterfly hysteresis loop observed in time-resolved magnetization measurements has been interpreted in terms 
of a combination of a Landau-Zener-Stiickelberg transition at zero field and spin-phonon coupling Our results 

show that unless the field is applied in a special direction (x or y-direction in this case), the adiabatic magnetization 
process is no longer symmetric with respect to the field. The dependence on the direction of the field should lead 
to observable changes in the hysteresis loops. So far, only weak directional dependence has been reported in experi- 
ments 44| . Therefore it seems that it is necessary to explore other mechanisms that yield energy level repulsions such 
as hyperfine interactions j34[. 
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Appendix A: Diagonalization of model (0 

Here we collect some analytical results of the solution of the eigenvalue problem of Hamiltonian Q). We only 
consider the case of a magnetic field that is parallel to the z-axis (h = (0,0, h z )). For a DM vector satisfying the 
conditions J2J) and (0 the eight eigenvalues are given by 



Eqa = 



-s/3D z + 4h z ±\9D 



9D V 



{V3D Z — 3 J + 2h z 



E- 



V3D Z - 4h z ± J9D X 2 + 9D y 2 + {V3D Z - 3 J — 2h z 



2.3 



E. 



V3D Z ± \ZD X 2 + 3D y 2 + {V3D Z + 3 J - 2h z ) 2 



4,5 



^6,7 — 



V3L> 2 ± a/3^ 2 + 3D y 2 + (V3D Z + 3 J + 2h z ) 2 



(10) 



where D x — Df 2 , D y = Dj' 2 , and D z = Df 2 . Substituting the values of all model parameters, we recover the results 
obtained by numerical diagonalization. For h z = there are four pairs of two- fold degenerate levels. Although it is 
possible to find analytical expressions for the case that the magnetic field is parallel to the :r-axis, the expressions 
themselves are rather long and not very illuminating. Therefore they are not given here. 
For a magnetic field parallel to the z-axis, a straightforward calculation shows that 



W|3/2,3/2) 
TL 2 |3/2, 3/2) 

H |3/2,l/2) 
H 2 |3/2, 1/2) 



- 3 ( J + 2 ^ |3/2,3/2) + ^ 



iD v 



9{(J + 2h z ) 2 + D 2 x + D 2 y ] 0/ox 3(D x + iD y )(V3D z +4h z ) i 

-^,3/2,1/2)- ^■ 4 +iD »» | S) , 
(3J + 2h z ) 2 + 3(D 2 X + D 2 y 
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|3/2, 1/2) _3(^ + i^, |a)| 



where \a) denotes the state \a) with all spins reversed and 



2V3 



(i-iV3)\ |TT> + (i + *V3)| UT)-2| TTI) 



(ii) 



(12) 



The expressions for H |3/2, -3/2), U 2 |3/2, -3/2), H |3/2, -1/2), and TL 2 |3/2, -1/2) are obtained from Eqs. (TTJ and 
(fT^fl by changing the sign of h z and D y and replacing |o) by 



\b) = 



2\/3 



(1 - iV3)| Til) + (1 + *V3)| ITI) - 2| lit) 



(13) 



Note that (a\b) = 0. From Eq. I|ll|) it follows that for the external field parallel to the z-axis, model Q does not 
allow transitions from the state with all spins up (down) to the state with two spins down (up). Therefore, if initially 
the system is in the state with all spins down, adiabatically sweeping the field from a large negative value to a large 
positive value will not yield the final state with all spins up. 



Appendix B: Numerical methods 



A theoretical description of quantum dynamical phenomena in the Mni 2 and V15 nanomagnets requires detailed 
knowledge of their energy-level schemes. Disregarding the fascinating physics of the nanomagnets, the calculation 



13 



of the eigenvalues of their model Hamiltonians is a challenging problem in its own right. Firstly, the (adiabatic) 
quantum dynamics of these systems is mainly determined by the (tiny) level repulsions. Therefore the calculation of 
the energy levels of these many-spin Hamiltonians has to be very accurate in order to bridge the energy scales involved 
(e.g. from 500K to sw 10 _9 K). Secondly, the level repulsions originate from the DM interactions that mix states with 
different magnetization. In principle, this prevents the use of the magnetization as a vehicle to block-diagonalize the 
Hamiltonian and effectively reduce the size of the matrices that have to be diagonalizcd. If a level repulsion involves 
states of significantly different magnetization (e.g. M z = — 10 and M z = 10) a perturbative calculation of the level 
splitting would require going to rather high order (at least 20), a cumbersome procedure. Therefore it is of interest 
to explore alternative routes to direct but accurate, brute-force diagonalization of the full model Hamiltonian. 

As a non-trivial set of reference data, we used the eigenvalues obtained by full diagonalization (using standard 
LAPACK algorithms) of the 10000 x 10000 matrix representing model JSJ 42] . For one set of model parameters, such 
a calculation takes about 2 hours of CPU time on an Athlon 1.8 GHz/1. 5Gb system. Clearly this is too slow if we 
want to compute the energy-level diagram, in particular if we want to estimate the structure of the level splittings. 
At the resonant fields we need the eigenvalues for many values of h. Furthermore, in the case of Vis this calculation 
would take about 30 times longer and require about 15 Gb of memory which, for present-day computers, is too much 
to be of practical use. 

We have tested different standard algorithms to compute the low-lying eigenvalues of large matrices. The standard 
Lanczos method (including its conjugate gradient version) as well as the power method [53l l54j| either converge 
too slowly, lack the accuracy to resolve the (nearly)-degenerate eigenvalues, and sometimes even completely fail to 
correctly reproduce the low-lying part of the spectrum. This is not a surprise: By construction these methods work 
well if the ground state is not degenerate and there is little guarantee that they will work if there are (nearly)- 
degenerate eigenvalues |53l l54| . In particular, the Lanczos procedure suffers from numerical instabilities due to the 
loss of orthogonalization of the Lanczos vectors [53l l54| . It seems that model Hamiltonians for the nanoscale magnets 
provide a class of (complex Hermitian) eigenvalue problems that are hard to solve. 

Extensive tests lead us to the conclusion that only the Lanczos method with full orthogonalization (LFO) [E3,03| an d 
the Chebyshev-polynomial-based projector method (CPP) discussed below can solve these rather large and difficult 
eigenvalue problems with sufficient accuracy. The former is significantly faster than the latter but using both gives 
extra confidence in the results. 



Lanczos method with full orthogonalization 

In the LFO, each time a new Lanczos vector is generated we exp licitly orthogonalize (to working precision) this 
vector to all, not just to the two previous, Lanczos vectors [53l l54| . With some minor modifications to restart the 
procedure when the Lanczos iteration terminates prematurely, after n steps this procedure transforms the n x n 
matrix H into a tri-diagonal matrix that is comparable in accuracy to the one obtained through Householder tri- 
diagonalization but offers no advantages |54j |. In our case we are only interested in a few low-lying eigenstates of H. 
Thus we can exploit the fact that projection onto the (numerically exact) subspace of dimension k (k <C n) built by 
the Lanczos vectors will yield increasingly accurate estimates of the smallest (largest) eigenvalues and corresponding 
eigenvectors as k increases. 

In practice, to compute the M lowest energy levels, the LFO procedure is carried out as follows. 

• Perform a Lanczos step according to the standard procedure. 

• Use the modified Gramm-Schmidt procedure to orthogonalize the new Lanczos vector with respect to all previous 

ones mm. 

• Compute the matrix elements of the tri-diagonal matrix. 

• At regular intervals, diagonalize the tri-diagonal matrix, compute the approximate eigenvectors ipi, fj,i = 
(ifii\H\ipi) and Af = (ifi\(H — /j,i) 2 \(pi) for i = 1,...,M, and check if all are smaller than a specified 
threshold. If so, terminate the procedure (the exact eigenvalue Ei closest to ^ satisfies /U, — A < E- L < /a; + A,-). 
If not, continue generating new Lanczos vectors, etc. 



Chebyshev polynomial projector method 



As an alternative to the LFO, we have used a power method [53, |54( based on the matrix exponential e tH 
Writing the random vector ^(0) in terms of the (unknown) eigenvectors {4>i} of H, we find 
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V(t) = e- tE ° [0 o (^o|*(O))+e-*( Bl - f! «V 1 (0 1 |vl/(O))+ e - t ^- Bo )0 2 (0 2 |*(O)) + ...] , (14) 

showing linii_ (00 $(t)/||$(t)|| oc <fro if (0o|^(O)) ^ 0. In this naive matrix-exponential version of the power method, 
convergence to the lowest eigenstate is exponential in t if E\ > Eq . 

The case of degenerate (Eq = E\ = ...) or very close (Eq w E\ w ...) eigenvalues can be solved rather easily by 
applying the projector to a subspace instead of a single vector, in combination with diagonalization of e tH within this 
subspace |55). First we fix the dimension k of the subspace by taking k equal or larger than the desired number of 
distinct eigenvalues. The projection parameter t should be as large as possible but nevertheless sufficiently small so 
that at least the first k terms survive one projection step. Then we generate a set of random initial vectors ^i(O) for 
i = 1, . . . , k and set the projection count n to zero. We compute the k lowest eigenstates by the following algorithm |55j 

• Perform a projection step + l)t) = e~ tHi l>i(nt) for i = 1, . . . , k. 

• Compute the k x k matrices. A = (*j((n + l)t)\e tH \^i((n + l)t)) = (*i((n + l)t)\^i(nt)) and B = (#i((n + 
l)i) | (rx + l)i)). Note that A is Hermitian and -B is positive definite. 

• Determine the unitary transformation U that solves the k x k generalized eigenvalue problem Ax = XBx. Recall 
that k is small. 

• Compute + l)t) = £* =1 U l ^ ] {{n + l)i) for i = 1, . . . , k. 

• Set *i((n + l)t) = + for i = 1, . . . , k. 

• Compute m = l)t)|fl"|*<((n + 1)*)) and check if Af = l)t)| (if - /ii) 2 |*i((n + l)t)) is smaller 
than a specified threshold for i = 1, . . . , k. If yes, terminate the calculation. If no, increase n by one and repeat 
the procedure. 

We calculate e~ tH '& by using the Chebyshev polynomial expansion method j5|| IH, H^, |(3(J • First we compute 
an upper bound R of the spectral radius of H (i.e., ||i?|| < R) by repeatedly using the triangle inequality [6(|. From 
this point on we use the "normalized" matrix H = (2H/R — l)/2. The eigenvalues of the Hermitian matrix H are 
real and lie in the interval [—1,1] |H3. IH^. Expanding the initial value \&(0) in the (unknown) eigenvectors 4>j of H 
(or H) we find 

*(t) = e -^( ) = e »5*(0) = X;e«^^<^|*(0)), (15) 

i 

where z = —tR. We find the Chebyshev polynomial expansion of by computing the Fourier coefficients of the 

function e zcos [6l|. Alternatively, since —l<Ej < 1, we can use the expansion e zEj — Io(z) + 2 Y^m=i Im(z)T m (Ej) 
where I m (z) is the modified Bessel function of integer order m |6j to write Eq. (|T3|) as 



*(*) = 



I {z)I + 2j2l m (z)T m (H) 



*(0). (16) 



Here, I is the identity matrix and T m (H) is the matrix- valued Chebyshev polynomial defined by the recursion relations 

T (ir)tf(0) = ¥(0), Ti(fl)*(0)=5*(0), (17) 

and 

T„ i+1 (#)*(0) = 2HT m (H)V(0) - T ro _ 1 (5)*(0) , (18) 

for m > 1. In practice we will sum only contributions with m < M where M is chosen such that for all 
m > M, \I m (z) / Iq(z)\ is zero to machine precision. Then it is not difficult to show that \\e~ /Iq(z) — I — 
^Em=i[^"i( z ) / ' Io(z)]T m (H)\\ is zero to machine precision too (instead of e~ fH we can equally well use e~ tH / Io(z) as 
the projector). 

Using the downward recursion relation of the modified Bessel functions, we can compute K Bessel functions to 
machine precision using of the order of K arithmetic operations j^, |g3 ■ A calculation of the first 20000 modified 
Bessel functions takes less than 1 second on a Pentium III 600 MHz mobile processor, using 14-15 digit arithmetic. 
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Hence this part of a calculation is a negligible fraction of the total computational work for solving the eigenvalue 
problem. Performing one projection step with e~ tH amounts to repeatedly using recursion (|18|) to obtain T m (B)^ (0) 
for k = 2, ... , M, multiply the elements of this vector by I m {z) and add all contributions. 
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